home
***
CD-ROM
|
disk
|
FTP
|
other
***
search
/
Pascal Super Library
/
Pascal Super Library (CW International)(1997).bin
/
MATH
/
NRPAS13
/
MEMCOF.DEM
< prev
next >
Wrap
Text File
|
1991-04-29
|
681b
|
32 lines
PROGRAM d12r9(input,output,dfile);
(* driver for routine MEMCOF *)
CONST
n=1000;
m=10;
TYPE
glnarray = ARRAY [1..n] OF real;
glmarray = ARRAY [1..m] OF real;
VAR
i : integer;
pm : real;
cof,wkm : glmarray;
data,wk1,wk2 : glnarray;
dfile : text;
(*$I MODFILE.PAS *)
(*$I MEMCOF.PAS *)
BEGIN
glopen(dfile,'spctrl.dat');
FOR i := 1 to n DO BEGIN
read(dfile,data[i])
END;
close(dfile);
memcof(data,n,m,pm,cof,wk1,wk2,wkm);
writeln('Coefficients for spectral estimation of SPCTRL.DAT');
writeln;
FOR i := 1 to m DO writeln('a[',i:2,'] =',cof[i]:12:6);
writeln;
writeln('b0 =',pm:12:6)
END.